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Abstract. We present results from a calculation of the QCD equation of state with 
two light (up, down) and one heavier (strange) quark mass performed on lattices with 

^ three different values of the lattice cut-off. We show that also on the finest lattice 

, analyzed by us observables sensitive to deconfinement and chiral symmetry restoration, 

I I respectively, vary most rapidly in the same temperature regime. 

1. Introduction 

OO , The equation of state (EoS) of strongly interacting elementary particles is one of the most 
fundamental non-perturbative quantities that numerical studies of lattice regularized 
QCD will be able to provide as input to the hydrodynamic modeling of the expansion of 
^ ! dense matter formed in heavy ion collisions. At least for vanishing chemical potential, 
OO ! which is appropriate for the conditions met in experiments at RHIC and LHC, lattice 
calculations of the EoS [H [21 [3] as well as the transition temperature [H O [6] can now 
be performed with an almost realistic quark mass spectrum. Calculations at different 
^ i values of the lattice cut-off allow for a systematic analysis of discretization errors and 
5^ I will soon lead to a controlled extrapolation of the EoS with physical quark masses to 
the continuum limit. 

We present here results from calculation of the EoS on lattices with temporal 
extent A^,- = 4 and 6 [3], which correspond to two different sets of the lattice cut-off, 
aT = 1/Nt-. We furthermore report on preliminary results for the EoS obtained on even 
finer lattices with temporal extent Nr = 8 [7j . We also show results for the temperature 
dependence of the strange quark number susceptibility and the chiral condensate. These 
observables are sensitive to different non-perturbative aspects of the QCD transition. 
Their temperature dependence suggests that the onset of deconfinement as well as 
the gradual melting of the chiral condensate are correlated and happen in the same 
temperature regime. 

All results presented here are obtained from numerical calculations in lattice 
regularized (2-1-1) -flavor QCD using staggered fermions. The calculations have been 
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Figure 1. The trace anomaly, (e — 3p)/T^ (left) calculated on lattices with temporal 
extent Nt = 6, 8 and the ratio of pressure and energy density as well as the velocity 
of sound obtained in calculations with the p4fat3 action on Nr — 6 (short dashes) [3] 
and 8 (long dashes) [7J. The temperature scale has been obtained from an analysis of 
the heavy quark potential from which the Sommer scale parameter has been extracted. 
Its value has been fixed to tq — 0.469 fm. 

performed with physical values for the strange quark mass and two degenerate light 
quark masses that correspond to a pion mass of about 220 MeV. 

2. The QCD equation of state 

Systematic studies of the QCD EoS are currently performed with two different versions 
of improved discretization schemes for staggered fermions, the asqtad and p4fat3 actions. 
Both actions are constructed such that they remove in the high temperature limit 0(a^) 
lattice discretization errors in bulk thermodynamic observables and reduce explicit 
flavor symmetry breaking effects through the introduction of so-called fat links. In 
the construction of these actions different strategies have been followed to deal with 
these lattice artefacts. Both actions have been used for some time to study the QCD 
equation of state on lattices with temporal extent iV^ = 4 and 6 [2], |3]. In a joint effort 
the hotQCD collaboration currently extends the studies of the EoS with these actions 
to lattices with temporal extent A^^ = 8 [8l|9]. Like in the earlier calculations performed 
with the p4fat3 action large spatial lattices are used {N„ = iN^-) to get close to the 
thermodynamic limit. 

In Fig. [1] (left) we show results for the trace anomaly, O'^'' = e — 3p in units of T^. 
The figure shows results from calculations performed on lattices with temporal extent 
Nr = Q and 8 using the asqtad as well as the p4fat3 actions. It can be seen that both 
discretization schemes lead to quite good agreement in a wide range of temperatures; 
although a closer inspection shows still a cut-off dependence of the results. They lead 
to a reduction of the peak height in Q^^^/T^, which is located at T ~ 200 MeV, and 
lead to a shift of the rapidly rising part of Q^^^ /T^ in the transition region to smaller 
values of the temperature. 

Cut-off effects as well as differences arising from both discretization schemes seem 
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to be largest in the vicinity of the maximum of (e — 3p)/T^. This trend has already 
been observed when comparing results obtained on lattices with temporal extent A^,- = 4 
and 6 [3]. The cut-off dependence observed in Qt^t^/T^ carries over to the calculation of 
pressure and energy density; the former is obtained by integrating over Qt^t^/T^ and the 
energy density is then obtained by combining results for p/T^ and (e — 3p)/T^. This is 
apparent in Fig. [1] (right) where we show the ratio p/e obtained with the p4fat3 action 
for three different values of the cut-off. Cut-off effects are still visible in the vicinity of 
the 'softest point' of the EoS, which is related to the peak position of (e — 3p)/T^. In 
the entire range of energy densities relevant for the expansion of dense matter created 
at RHIC, e < 10 GeV/fm^, the ratio p/e deviates significantly from the conformal, 
ideal gas value p/e = 1/3. This also is reflected in the behavior of the velocity of 
sound, Cg = dp/de, which is shown in Fig. [T] (right) by dashed lines. It starts deviating 
significantly from the ideal gas value below e ~ 10 GeV/fm^ and reaches a value of 
about 0.1 in the transition region at energy densities e ~ 1 GeV/fm^. 

3. Deconfinement and chiral symmetry restoration 

The relation between deconfinement and chiral symmetry restoration in QCD has been 
discussed since a long time. Although both phenomena seem to be related to physics on 
different length scales lattice calculations seem to suggest that both phenomena happen 
at approximately the same temperature even at finite, non-zero values of the quark 
masses when none of the symmetries related to confinement {Z{3) center symmetry) or 
chiral symmetry breaking {SUL{nf) x SUji{nf)) are realized exactly. This aspect of the 
QCD transition received renewed interest recently; as the QCD transition is a crossover 
rather than a genuine phase transition different observables, that are sensitive to 
different aspects of the QCD transition, might lead to different transition temperatures. 
In particular, it has been suggested that observables sensitive to deconfinement may 
lead to a higher transition temperature than observables sensitive to chiral symmetry 
restoration [6]. 

With decreasing quark mass the Polyakov loop looses its role as a genuine order 
parameter for deconfinement. It is non-zero at all temperatures. Nonetheless it varies 
rapidly in the transition region indicating that quark free energies are screened more 
effectively in the high temperature phase of QCD [3]. Another observable, sensitive to 
deconfinement is the quark number susceptibility. 



where q = I, s for the light and strange quark sector, respectively. These susceptibilities 
are sensitive to the degrees of freedom that carry a net number of light or strange quarks, 
respectively. They too change rapidly in the transition region as the carrier of, e.g. 
strangeness, are rather heavy strange hadrons at low temperature and single quarks at 
high temperature. A rapid rise in the susceptibilities thus indicates 'deconfinement'. 
In the limit of massless quarks the chiral condensate, {'ip'ip)q = TV^^dln Z/drUq, is 



T2 VT^ difXg/T) 



2 ' 



(1) 
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Figure 2. The strange quark number susceptibility (left) and the subtracted chiral 
condensate normalized to the corresponding zero temperature value (right). The band 
corresponds to a temperature interval 185 MeV < T < 195 MeV. Results shown for 
Nr = 8 are preliminary data of the hotQCD collaboration. 

an order parameter for chiral symmetry restoration. At non-zero values of the quark 
mass the condensate receives additive as well as multiplicative renormalization which 
have to be eliminated to allow for a sensible continuum limit. An appropriate observable 
is obtain through subtraction of a fraction of the strange quark condensate from the light 
quark condensate. The difference, taken at finite temperature, can then be normalized 
with the corresponding zero temperature difference. 
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This observable, in addition, has a sensible chiral limit and gives, in this limit, an order 
parameter that is unity at low temperature and will vanish at Tc. 

In Fig. [2] we show the strange quark number susceptibility (left) as well as Ai^s{T) 
(right). It is apparent from these figures that observables sensitive to deconfinement and 
chiral symmetry restoration, respectively, vary most rapidly in identical temperature 
intervals. This interval coincides with the region of most rapid rise in the trace anomaly 
as well as energy and entropy density. 

This manuscript has been authored under contract number DE-AC02-98CH1-886 
with the U.S. Department of Energy. 
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